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Abstract. The standard model for radio supernovae considers that the observed synchrotron radio emission arises 
from the high-energy shell that results from the strong interaction between the expanding supernova ejecta and 
the circumstellar medium. This emission is considered to be only partially absorbed by ionized thermal electrons 
in the circumstellar wind of the progenitor star. Based on a study of the radio light curves of the type II supernova 
SN1993J, we present evidence of synchrotron self-absorption. Our modeling of the radio light curves requires a 
large initial magnetic field, of about 30 Gauss, and the existence of an (initially) highly-relativistic population of 
electrons. We also show that while at early epochs the dominant absorption mechanism is external absorption 
by thermal electrons, at late epochs and long wavelengths the dominant absorption mechanism is synchrotron 
self-absorption. Consequently, the spectral turnover takes place at much shorter wavelengths than expected in the 
standard model, and at long wavelengths (> 90cm at current epochs) the flux predictions depart substantially 
from those of the standard model. 
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1. Introduction 

Radio emission from supernovae has been previously in- 
terpreted in terms of an optically-thin shell, whose emis- 
sion is partially suppressed by external free-free absorp- 
tion (e.g., Chevaher |l982a| ; Weiler et al. |1989| , |1990D . In 
this model - known as the standard interaction model 



(Chevalier 1982a), hereafter SIM - a strong interaction be- 
tween the expanding supernova ejecta (pej cx r~") and the 
circumstellar medium (pcs oc r~^) is expected. In the SIM, 
this interaction causes the formation of a self-similarly 
expanding {Rsh oc t™;m — {n ~ 3)/(n — s)) shell-like 
structure (Chevalier [l982b| , Nadyozhin |l98|) , fr om which 
the observed synchrotron radio emission arises. In addi- 



tion, it is assumed (Chevalier 1982a) that (i) both mag- 
netic energy density, ee (x B^, and relativistic energy den- 
sity, Crei cx Nq E^~p , cvolvc with time as the post-shock 
thermal energy density (oc Psh^sh)' (^0 synchrotron 
emission is optically thin; and (iii) the external absorbing 
medium has a power-law dependence s = 2. 
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The discovery of SN 1993J (Garcia |1993|) on March 28, 
1993 in the galaxy M81, brought an unprecedented oppor- 
tunity to test the SIM for supernovae. The early discovery 
of hy drogen lines in its spectra (Fillipenko & Matheson 
1993 ) classified this supernova as Type II. Marcaide and 
coauthors showed the radio emission to arise from a shell- 



like structure (Marcaide et al. 1995a ) that expands self- 
similarly (Marcaide et al. 1995b ), thus g iving strong sup- 
port to the SIM. Marcaide et al. ( 1997 ) also reported on 
the deceleration of the expansion (m = 0.86 ± 0.02), and 
concluded that there is strong evidence for the circumstel- 
lar medium density profile to be of the form pcs cx r~*, 
with s = 1.66, in good agreement with previous model- 
ing of radio observations (s = 1.5; Van Dyk et al. 1994, 
Lundqvist 1994) and with the modeling of both radio 
and X-ray observations, which require 1.5 ^ s 1.7 
(Fransson et al. 1996; hereafter FLC96). Attempts to fit 
the radio light curves of SN 1993 J using the SIM have 
been unsuccessful so far, and have led several authors (Van 
Dyk et al. |1994[ Lundqvist |1994 FLC96, Marcaide et al. 



1997) to consider alternative power-law profiles for the cir- 
cumstellar medium around SN 1993J, or to include syn- 
chrotron self-absorption as an additional absorbing mech- 



anism (Fransson & Bjornsson 1998, hereafter FB98). 
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2. Synchrotron radio emission from SN 1993 J 

In an attempt to gain insight into the physical conditions 
of SN 1993 J, and also motivated by our VLBI monitor- 
ing of this unique radio supernova, we have written a nu- 



merical code (Perez- Torres 1999) that computes the syn- 
chrotron radio emission from an expanding Type II super- 
nova. Although our code will be fully described elsewhere, 
we briefly outline it here. We solve the synchrotron radia- 
tive transfer equation throughout the shell (including syn- 
chrotron emission and synchrotron self-absorption), and 
take into account both synchrotron and expansion losses. 
We assume a (modified) SIM scenario, in which the su- 
pernova synchrotron radio emission comes from the inter- 
action shell that forms shortly after shock breakout. (The 
SIM scenario is modified in the sense that we do not as- 
sume the synchrotron emission to be optically thin.) We 
assume the emission to be due to the existence of rela- 
tivistic electrons, N(E^ r), orbiting in a random magnetic 
field, B{r). We assume an initially injected relativistic 
electron distribution with a power-law energy dependence, 
N{E, r) — NQ{r) E^p, and also assume a continuous reac- 
celeration of those electrons. 

In modeling the radio emission from SN 1993J, we take 
into account all available observational data: (i) the dis- 
tance to SN1993J, D = 3.6 Mpc (Freedman et al. |1994D; 



(ii) the supernova explosion date, icxp = March 28.0, 1993 



(Nom oto et al. |1993| , Podsiadlowski et al. |1993 , Bartunov 
et al. 1994); (iii) the supernova deceleration parameter, 
0.86 (Marcaide et al. |1997| ); (iv) the width of the 



radio emitting shell, AR = 0.3 i?shock (Marcaide et al. 



1995b, 1997); (v) the index of the circumstellar density 
material, s = 1.66 (Marcaide et al. 19971 ). Finally, we 



assume (vi) the same external electron temperature pro- 
file as in FB98: Te{r) = max [T15 (lO^^ cm/r) , 2 x 10^ K] , 
where Tis is the electron temperature in the circum- 
stellar medium at r = 10"'^^ cm, which in our model 
is equal to 1.7 x 10® K. The thermal electron density 
of the circumstellar medium (assumed to be fully ion- 
ized with abundances X = 0.73,1^ = 0.23), is given 
by Ties = 1.27 X 10^ Mw (r/6.35 X lO^-^cm)"" cm~3^ 
where we have introduced the wind parameter Mw = 
(^M/lQ-^Moyr-i) (w„/10kms-i)~^ gcm-i, and M 
and Vw are the mass-loss rate and velocity, respectively, 
of the presupernova wind. As in the SIM, we assume that 
both B and Nq evolve with time as the post-shock ther- 
mal energy density, PshVgh- With m = 0.86 and s = 1.66 
from our VLBI observations, we then obtain the fol- 
lowing evolution laws for magnetic field and relativistic 
particles, respectively: B cx 7.i-''/2-i/m _ 7.-0.99^ ^ 

^2(l-s/2-l/m) ^ ^-1.99^ ^-^^^ jg^ ^ ^ ^-1^ g^^^ ^ ^^2 _ 

These scalings correspond to those in model 4 of Chevalier 
(|1996D , which give the slowest evolution of B and Nq, 
and therefore of the radio emission from SN1993J. The 
other three models consider a dependence of _B oc t^^, or 
No (X or both. 




Time since explosion (days) 



Fig. 1. Radio hght curves for SN 1993J at A1.3 cm (filled 
circles), A2 cm (open circles), A3. 6 cm (filled triangles), 
A6 cm (open triangles), and A21 cm (filled squares), from 
measurements made using the VLA (Van Dyk et al. 1994 , 
and Van Dyk, private communicatio n) an d, at 2 cm, also 
the Ryle Telescope (Pooley & Green 1993| ). The lines rep- 
resent our best fit model as described in the text. 



We model the available radio data using the following 
parameters; (1) the spectral index of the injected elec- 
tron distribution, p; (2) the initial magnetic field, i?(ro); 
(3) the initial value of the injected function of relativis- 
tic electrons, A'^o(''o); and (4) a low-energy cut-off for 
the relativistic electrons, 7min(''o) = Emin{ro)/meC^- We 
fit the above parameters at the reference radius tq — 
6.35 X 10^"* cm, corresponding to the outer shell radius 
at our first epoch {to = 2.34 days). We find that the 
multi-wavelength radio data for SN 1993J (Figure ^ are 
best reproduced by a model with p ~ 3, B{ro) = 27 G, 
No{ro) = 6.7 X IQ-^ergP"! cm'^, and 7min(ro) « 90. 

We estimate the uncertainties in the parameters of our 
model by modifying each of the best-fit parameters (one at 
a time) , and comparing via a least-squares fit the resultant 
emission with that predicted by our best-fit model. These 
calculations show the extreme sensitivity of our results 
to p, which needs to be kept very close to 3 (within a 
small few per cent) to obtain any meaningful fit to the 
radio light curves. The other three parameters are less 
well constrained, and a 20% variation of each of their 
best-fit values result in reasonable fits. In addition, we 
have considered the existence of couplings between some 
of the parameters, especially No and B. Further testing 
shows that factor of two variations in each of the three 
poorly-constrained parameters give very poor fits to the 
observations. Therefore, we adopt factor of two modeling 
errors for each of those three parameters. 

2.1. Spectral evolution of SN 1993 J: evidence for 
synchrotron self-absorption 

Our simulated spectra for SN 1993J agree well with the ob- 
served spectra (Fig. ||). The supernova spectrum changes 
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Fig. 2. Fits to spectra at 11, 237, and 1397 days after ex- 
plosion, and predicted spectrum at 2700 days, considering 
external free-free absorption (dashed lines), synchrotron 
self-absorption (dotted lines), and both external free- free 
absorption and synchrotron self-absorption (solid lines). 
The filled circles represent the available observational data 
for each epoch. At early epochs, the radio spectrum of 
SN 1993J (upper left plot) is well explained by the solely 
effect of free- free absorption from the external medium. At 
late epochs (lower left plot), however, both free- free ab- 
sorption and synchrotron self-absorption predict a similar 
spectral behavior for the observed wavelength range, but a 
different one at wavelengths around and above A40 cm, at 
which observations are not available. The lower right-hand 
plot shows the predicted spectrum at current epochs. This 
prediction can be tested to discriminate between the (ex- 
ternal) free-free and (internal) self-absorption mechanisms 
described in the text. 



from a mainly inverted spectrum at early epochs to a more 
typical optically thin spectrum at late epochs and for an 
increasingly wide wavelength range. As shown in Fig. H, 
external free-free absorption is dominant at early times 
{t < 100 days), but synchrotron self-absorption becomes 
increasingly important with time for the long-wavelength 
radio emission. Thus, the standard interaction model is 
incomplete, at least for SN 1993J, in that it ignores syn- 
chrotron self-absorption, which according to our modeling 
is the most relevant absorption mechanism - as the super- 
nova ages - for the longer wavelengths. At current epochs, 
our model predicts that the flux at A90cm is ^^120 mJy for 
the synchrotron self-absorption dominated scenario, and 
~200 mJy for the free- free absorption dominated scenario. 
This prediction - which would represent a direct test to 
discern the relevant absorption mechanism for SN 1993J - 
can be checked using A ~90 cm observations of SN 1993J 
with a sensitive radio interferometer. 



3. Results and Discussion 

Our modeling results show that a combination of syn- 
chrotron self-absorption and external free-free absorption 
can reproduce well the evolution of the radio light curves 
of SN 1993J. Therefore, the ad hoc external inhomoge- 
neous medium ( "clumpy external absorption" ) previously 
invoked by Van Dyk et al. ( p_994 ) to model the radio light 
curves seems to be unnecessary. We have also considered 
the Razin-Tsytovich effect (the suppression of the emis- 
sion below a certain wavelength due to the existence of a 
plasma) as an alternative mechanism to free-free absorp- 
tion. However, the densities of thermal electrons required 
for the Razin-Tsytovich mechanism to be effective are un- 
realistically large. For example, at 11.5 days, we find that 
the number density of thermal electrons should be larger 
than ~ 1.3 X lO^cm^'^ for the Razin-Tsytovich effect to 
be important. This value is about an order of magnitude 
larger than needed to correctly model the X-ray lumi- 
nosity at those epochs (FLC96). The same holds at late 
epochs, and we therefore rule out the Razin-Tsytovich ef- 
fect as a plausible alternative to external free-free absorp- 
tion. We point out here that the index of the circumstel- 
lar medium density profile, s = 1.66, was obtained by 
Marcaide et al. (1997) under the assumption that exter- 
nal free-free absorption was the only absorption mecha- 
nism acting in SN1993J. Although there is increasing evi- 
dence that synchrotron self- absorption is also relevant for 
SN 1993J (Chevalier [l99^ , FB98, this work ), we note that 
at the early epochs used by Marcaide et al. (1997) free- free 
absorption is the dominant mechanism, while synchrotron 
self- absorption can be neglected (see Figure 2, epoch 11 
days). We have tried to fit the light curves using either 
an standard index for the circumstellar wind (s = 2), or 
a much flatter one (s = 1.5) which would still satisfy the 
range of 1.5 s 1.7 inferred from X-ray observations 
(FLC96). As Fig. ^ shows, an s = 2 wind causes a fall-off 
of the observed emission which is far too rapid. On the 
contrary, a wind with s = 1.5 makes the radioemission 
to fall-off too slowly. For the case s = 2, we also investi- 
gated different scalings for B and iVo; namely, w e use d 
the other three scalings mentioned by Chevalier( |l996| ) . 
However, none of those scalings gives good fits because 
all of them make both B and TVq -and therefore the ra- 
dio light curves- to decrease more rapidly than r^^ and 
r~^, respectively. This suggests that the index s should be 
close to our nominal value of s = 1.66. Since the decline 
rate of the radio light curves, i.e., their late time emission, 
is also affected by the parameter p, we have investigated 
models in which both s and p are allowed to vary. This 
search resulted in much poorer fits than our best one. We 
acknowledge, though, that the assumption of ad hoc scal- 
ings for both B and A^o could alleviate these poor fits, at 
the expense of more complicated models. We defer such 
detailed study to a further publication. 

We investigated as well the effects of using different in- 
dices for the expansion (to), but do not find a reasonable 
fit to the supernova emission using either an undeceler- 
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Fig. 3. Observed radio light curves for SN 1993J at 1.3, 
2, 3.6, and 6cm, and model predictions for three different 
indices of the presupernova circumstellar wind: s = 1.5 
(dotted lines), s = 1.66 (solid hnes), s = 2 (dashed hnes). 
See text for details. 



ated expansion of the supernova (m = 1), or a heavily 
decelerated expansion (m = 0.75). 

The unexpected large (initial) magnetic field of about 
30G implied by our modeling agrees with the large mag- 
netic field obtained by FB98, and seems to rule out any 
hypothesis based on a magnetic field due to a compression 
of the circumstellar wind field. Indeed, since it is very un- 
likely that the magnetic field energy density in the wind is 
larger than its kinetic energy density, i.e., B'^/Stt < pv'^/2, 
we should expect a magnetic field not in excess of i3 < 3.4 
mG at r ~ 3 X lO^^cm (if it arises from a compression of 
the magnetic field in the wind). This value of the magnetic 
field is in agreement with the values obtained from polar- 
ization observations of OH masers in supergiants (e.g., 
Cohen et al. |1987| , Nedoluha & Bowers [19 9 2] ), which give 
magnetic field estimates at distances of ~ 3 x lO^^cm of 
1 — 10 mG. However, we find that the magnetic field is 
actually about 0.58 G, i.e., 170 times larger. From a strong 
supernova shock (factor of 4 compression) , one would ex- 
pect an increase of the field value of ~ 4^^^ for a magnetic 
field that scales as B (x r^^ oc nl^^ , clearly insufhcient 
by almost two orders of magnitude. Therefore, the large 
magnetic field argues strongly in favor of turbulent am- 
plification in the she ll (e.g .. Gull 1973, Chevalier 1982a, 



Chevalier & Blondin |1995[) . Furthermore, the assumption 
that the magnetic energy density evolves as the post-shock 
thermal energy results in a dependence of the magnetic 
field with radius of r~^. A component of the field be- 
having like B oc r"^, as suggested from magnetic flux 



conservation, does not give a reasonably good fit to the 
observations. 

We find that the magnetic energy density in the shell 
is just a (small) fraction of the post-shock thermal en- 
ergy density, es/^sh ~ 0.024, but nonetheless much larger 
than the relativistic electron energy density, eid/esh ~ 
2 X 10~^ — 10~^. This points to a strong deviation from 
equipartition. We note, however, that the ions (whose 
emission we neglect as it is a factor nip/me « 2000 less 
than that due to the electrons) might have a large con- 
tribution to the total energy density. We also find that 
the relativistic particle number density is a very small 
fraction of the particle thermal density. In fact, the ra- 
tio nrci/n,h = (iVo-Bni7n/2"sh) « (1 - 8) x 10"^ at any 
given epoch. 

Special attention deserves the low-energy cutoff at 
early epochs that we consider to account for the flatten- 
ing of the radio light curves at early epochs. We recall 
that our model assumes an initial injection of (highly) 
relativistic particles N{E,r) = No{r) E^p , which applies 
from Emin up to Emax- The initial supernova explosion 
is the likely seed for such a population of relativistic elec- 
trons, which must be continuously reaccelerated. We show 
in Fig. H the light curves that we obtain when we do 
not consider a low-energy cutoff in the relativistic parti- 
cle distribution, i.e., 7min = 1- (The other parameters are 
fixed to their nominal values p — 5, B{ro) — 27 G, and 
iVo(ro) = 6.7 X lO-'^ergP-icm-^.) We note that although 
the optically thin part of the radio light curves remains un- 
affected, the optically thick part is heavily modified by the 
lack of a low-energy cutoff. In particular, the radio emis- 
sion rises more rapidly. More importantly, the observed 
flattening of the light curves at the shortest wavelengths 
is not reproduced by such a model. Notice also that at long 
wavelengths (see the light curve at A21cm in Fig. ^ the 
radio emission is strongly absorbed at early times. This 
is the result of having decreased the minimum energy of 
the electrons; since the number of relativistic particles is 
proportional to E^^^, we effectively increased the column 
density of the relativistic particles by a large factor. The 
effect is most easily seen at long wavelengths due to the de- 
pendence of the opacity with wavelength {ki^ oc A'p^**'/^). 
At this point, we also note that our model differs here from 
that presented by FB98, in that they assume a continu- 
ous injection of particles with 7 > 1. Bell (1978) showed 
that particles that are already relativistic can be accel- 
erated in the presence of a shock. Since it is usually as- 
sumed that the acceleration of particles in SNe is due to 
the Fermi process (e.g. Ball & Kirk |1995| for SN1986J), we 
show here that our finding of a low-energy cutoff is consis- 
tent with such a process acting in SN1993J. The minimum 
energy of electrons that can be efficiently accelerated is 
7min,c = (TTir, I rne) {v A I c iA (Eilck & Hughcs 1991; but see 
also Ellison et al. 200C| , who argue that nonlinear effects 
in shocks can lead to acceleration of particles from a ther- 
mal pool), where nip and vrip are the protron and electron 
mass, respectively, Vf^ is the Alfven velocity, c is the speed 
of light, and /i is the cosine of the pitch-angle. With the 
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Fig. 4. Radio light curves for SN 1993J. The data is as 
in Fig. 1^. The lines correspond to our best-fit model as 
described in the text, except that we did not assume a low- 
energy cut-off. Such model fails to reproduce the flattening 
of the short-wavelength emission seen at early epochs. See 
text for details. 



physical conditions deduced for SN 1993J, va varies from 
« 3200 km s^^ at t=2.3 days down to about 1200 km s^^ 
at t=3000 days. This gives 7i„in > 20 and 7„iin > 7 at 2.3 
and 3000 days, respectively, which is consistent with our 
results. 

It is plausible that Coulomb losses also play a role in 
the shape of the SN 1993J light curves at high frequen- 
cies and early epochs (see FB98, which use power-laws 
for the circumstellar wind density and the deceleration 
of the expansion of SN 1993J that are different fr om th e 
experimental values obtained by Marcaide et al. 1997 ). 
However, the neglect of Coulomb losses in our modeling 
is self-consistent, and does not affect the behavior of the 
radio light curves either at early or at late epochs. Indeed, 
the Lorentz factor, 7, of the electrons responsible for the 
emission at a given wavelength is 7 « 84 (A i? sin 9) 
(cgs units), while the Lorentz factor below which Coulomb 
losses dominate over expansion losses is approximately 
7Coui « 180 (nsh/6.9 x lO^cm'^) (^/10days)^"™^ where 
we have neglected for simplicity the slowly-dependent 
term on the electron energy. From our modeling at t = 10 
days, we get B « 7.7G and 7min ~ 61, so the ob- 
served emission at A 1.3 cm comes mainly from electrons 
with 7 « 200, while Coulomb losses are important for 



7 ^ 7Coui ~ 180. Analogously, our modeling gives a mag- 
netic field of ^ 65 mC at current epochs (~2700 days), 
and 7niin « 12. We then obtain that the observed A90 
cm emission comes mainly from electrons with 7 ^120, 
while Coulomb losses are relevant only for electrons with 
7 ^ Icoui ~ 17, and thus do not affect the long wavelength 
radio emission from SN 1993J at late epochs. Therefore, 
we stress that our main result, namely the dominance of 
synchrotron self-absorption over external free- free absorp- 
tion at late epochs is a robust result, and is not affected by 



the physical details that cause the flattening of the spectra 
at early epochs. 

In summary, we have presented evidence of syn- 
chrotron self-absorption for the type II SN 1993J, based on 
a modeling of its radio light curves. We show that the radio 
light curves of SN 1993J are better characterized by syn- 
chrotron emission that is partially absorbed by a combina- 
tion of synchrotron self-absorption (inside the shell) and 
external free-free absorption, rather than only by external 
free-free absorption, thus confirming previous modeling of 
the SN 1993 J radio hght curves (FB98). This result sug- 
gests that the "clumpy external absorption" previously in- 
voked to model the radio light curves of SN 1993J might 
not be needed. We also find that the Razin-Tsytovich ef- 
fect can be ruled out as well, as an alternative mecha- 
nism to external free- free absorption. We predict that syn- 
chrotron self-absorption is the main absorbing mechanism 
for the long-wavelength radio emission as SN 1993J ages. 
This prediction can be tested unambiguously using cur- 
rent technology and observations of SN 1993J at A90 cm. 
Moreover, follow-up observations of SN 1993J at A > 90 
cm are also needed to further constrain the physical pa- 
rameters of this unique supernova. As pointed out earlier 
(Chevalier |l998| ), if SN 1993 J is not a bizarre supernova 
then synchrotron self-absorption is likely to play an impor- 
tant role in most, if not all, radio supernovae, especially 
as they age. 
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